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Abstract 



We consider the interaction of two skyrmions in the framework of the sudden approx- 
■ imation. The widely used product ansatz is investigated. Its failure in reproducing an 

attractive central potential is associated with terms that violate G-parity. We discuss the 
construction of alternative ansatze and identify a plausible solution to the problem. 

Keywords: Skyrme model, nucleon-nucleon interaction, chiral symmetry, effective La- 
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1 Introduction 

Nucleon-nucleon (NN) interactions are relatively simple at large distances and become rapidly 
more complex as one moves inward. In the best phenomeno logical models existing at present, 
that reproduce low-energy observables accurately, they are described by the consensual one- 
pion exchange potential (OPEP), supplemented by theoretical two-pion exchange potentials 
(TPEP) and parametrized at short distances |], |2|. The OPEP is responsible for a strong ten- 
sor component, which is mostly important in few-body systems, such as the deuteron. Two-pion 
exchange, on the other hand, gives rise to the central potential, that survives to all averages 
and is responsible for most properties of large systems and nuclear matter. Quantum Chro- 
modynamics (QCD) is the basic framework for the study of strong processes and should have, 
in principle, an important role in the description of nuclear forces. However, at present, the 
non-Abelian character of this theory prevents low energy calculations and one has to resort 
to effective theories, which must reflect the main features of QCD. Thus, in Nuclear Physics 
applications, besides the usual space-time invariances, one requires these theories to have ap- 
proximate chiral symmetry. The latter is usually restricted to the SU(2) x SU(2) sector, for 
most processes involve only the quarks u and d. This symmetry is explicitly broken by the 
small quark masses and, at the effective level, by the pion mass. 
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Chiral symmetry has no influence over the OPEP, but is crucial to the TPEP, which depends 
on an intermediate pion-nucleon (ttN) amplitude [I]. In the case of NN interactions, the 
importance of this symmetry was stressed already in the early seventies, by Brown and Durso 
[[| and by Chemtob, Durso and Riska Q, who used it to constrain the form of the TPEP. 
In that decade it also became popular to describe nuclear processes by means of the linear a 
model H, containing a fictitious particle called a that, to some extent, simulates the TPE. 
The elimination of this unobserved degree of freedom gave rise to non-linear theories, which 
underlie modern descriptions of the interaction. The first theoretical framework to incorporate 
non-linear chiral dynamics into the NN problem was proposed by Skyrme. This remarkable 
model for the nucleon, developed in the sixties || and revived in the eighties @], describes 
baryons as topological solitons, objects extended in space that rotate according to the laws of 
Quantum Mechanics. The quark condensate appears as an intrinsic feature, corresponding to 
a non-vanishing classical content of the vacuum, whose intensity is given by the pion decay 
constant f^. Skyrmions correspond to distortions of this condensate that carry topological 
charges. One then works with pion fields which are unusually strong, in the sense that their 
amplitudes may be comparable to f^. Thus, in spite of its well know limitations ||, the Skyrme 
model remains a unique laboratory for studying chiral symmetry in the non-perturbative regime. 

In the early nineties Weinberg restated the role of perturbative chiral symmetry in nuclear 
interactions |§ and motivated interest in the TPEP. Initially, several authors explored the 



pion-nucleon sector of non- linear Lagrangians [10 , but the corresponding potentials could not 



reproduce even the medium range attraction in the scalar channel. This happened because the 
TPEP is based on an intermediate ttN amplitude, that can only be well described with the help 
of other degrees of freedom pT| . Accordingly, in a later stage, agreement with empirical ttN 
information was enforced and descriptions could reproduce NN scattering data [|1^, [13|, |T4|| . 

In the case of perturbative calculations, the delta is by far the most important non-nucleonic 
degree of freedom and is largely responsible for the intermediate range scalar attraction. As 
the Skyrme model incorporates the delta from the very beginning, one expects that it should 
yield a good qualitative NN potential. However, it fails to do so. 

Skyrme himself considered NN interactions, already in the sixties, using the so called 
product ansatz (PA) |J. The basic idea underlying the PA is that solutions corresponding 
to baryon number B = 1 can be used as building blocks to construct approximate solutions 
with an arbitrary value of B. The great advantage of this approach is that the baryon number 
of the composite system is automatically equal to the number of individual B = 1 skyrmions, 
irrespectively of their relative positions. In the PA, the skyrmions that constitute a larger 
system are assumed to retain their shape all the time, what is known as sudden approximation. 
In this framework, the construction of the NN potential is rather simple and the fact that each 
nucleon has a profile function which falls off rapidly with the distance allows one to assume 
that, for medium and large distances, the B = 2 system is not considerably different from the 
superposition of two with B = 1. 

In the eighties, the product ansatz was used by Jackson et al. []T5| and Vinh Mau et 



al. |L6| to calculate the NN potential, who found out a fully repulsive central component, 
in disagreement with very well established phenomenology. This puzzle motivated several 
attempts to construct improved versions of those early works. Among them, one notes the 
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symmetrized product ansatz by Nyman and Riska flT7| |, which could produce an intermediate 
range scalar attraction. However, as pointed out by Sternheim and Kalbermann [|1|, there is 
a violation of baryon number conservation in this ansatz. Exact numerical calculations were 
also used, which allowed one to evaluate the reliability of the sudden approximation at short 



distances [19]. Lattice calculations, using a method developed by Manton and collaborators, 



gave rise to a torus-like baryon density, believed to correspond to the true B = 2 ground 
state and having almost twice the nucleon mass [pOfl . The scalar potential associated with 
this configuration does show some medium to long range attraction JH]]. However, it is worth 
recalling that lattice results depend on the definitions adopted for collective coordinates and, 
as the full treatment is rather difficult, one usually resorts to approximations [22, 23| . 

In this work we consider the scalar interaction between skyrmions, in order to explore 
the possibility of obtaining the central attraction at large distances by relaxing some of the 
constraints present in usual calculations. We employ the sudden approximation because it gives 
rise to a constructive interaction, in which undeformed nucleons are the main building blocks, 
as in perturbative calculations. Our presentation is divided as follows. In section |^ we study the 
asymptotic behaviour of the scalar potential in the standard product ansatz approximation, in 
order to understand why it does not yield attraction. In section |3| we discuss the construction 
of alternative solutions, which must be constrained to have the correct baryon number. Finally, 
in section § we analyze a possible solution to the problem, and present concluding remarks. 



2 Central potential 



The structure of the central potential has been studied recently, in the framework of chiral 

m 



perturbation theory 
form 



V c (t) = - 



. In momentum space, the leading contribution has the generic 
2 
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where f w is the pion decay constant, the are subthreshold coefficients |Tl[ and a(m; t) is the 
scalar form factor, that depends on both the momentum transferred t and the baryon mass m. 
This form factor is defined generically in terms of the symmetry breaking Lagrangian jC s b as 



(N(p')\-C sb \N(p)} 
In configuration space, eq. ([I]) becomes 



V c {d) 



ft (c 



00 



a(t) u(p') u(p) 



+ doi V 2 J cr(m; d) 



(2) 



(3) 



where d is the internucleon distance and a(m; d) is the Fourier transform of a(m; t). In order 
to allow this result to be compared with the corresponding one in the Skyrme model, we note 
that, in the large N c limit, the nucleon and the delta are degenerate and very heavy. In ref. 



23J, the heavy baryon limit of a(m; d) was considered and found to be 
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with x = m w d. The central potential is 
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and, at very large distances, it behaves as 
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The sign of the constant K is determined by the values of the subthreshold coefficients in the 
combination (aQ + 4m 2 r a f L ). In table [l] we display empirical values for the and it is possible 
to note that the correct sign of Vc comes mainly from a^, since a^o in isolation would give rise 
to a repulsive interaction. 





ref. 11] 


ref. |H 


ref. [H 


ref. |27| 


4o ("O 


-1.46 ±0.10 


-1.30 ±0.02 


-1.27 ±0.03 


-1.15 ±0.03 




1.14 ±0.02 


1.35 ±0.14 


1.27 ±0.03 


1.23 ±0.03 


K(MeV) 


21.6 


28.5 


26.5 


26.2 



Table 1: Empirical values for the subthreshold coefficients ago, Oqi and the constant K, which 
determines the intensity of the central potential. 

In order to study the central potential in the Skyrme model, we recall that the standard 
soliton Lagrangian density is written as || ||] 

£ = C a + U , (7) 

where 

f 2 f 2 
C a = ^-Trid^Ud^tf) Tr(C7 + U ] - 2) (8) 

corresponds to the non-linear a model and 

U = -^Tr[d»Urf,d v Urf]* (9) 

is the stabilizing term. In these expressions, e is a free parameter, called Skyrme constant, 
whereas the dynamical variable U is a 2 x 2 unitary matrix, given by 

U = e lT - izF = cosF±ir-7rsinF, (10) 

where r are the isospin Pauli matrices and F is the chiral angle, whose boundary conditions 
determine the baryon number of a particular configuration. The function F and the isospin 
direction 7r are related to the pion field 7r of the non-linear a model || by n = f n sin F tz. 
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In the B = 1 static solution is obtained using the condition tt = r, the so called 

hedgehog ansatz, with boundary conditions F(r = 0) = tt and F(r — > oo) = ||. The 
quantization of this baryon is achieved by rotating the static solution with the help of collective 
coordinates, as a rigid body. This procedure endows the skyrmion with spin and isospin and 
corresponds to multiplying the pion field by the rigid body rotation matrix D, 

TTi -> K = D a i TTi- (11) 

The matrices D satisfy the completeness relations D ai D a j = 5^, D ai Dpi = 5 a/ 3 and, in the case 
of nucleons, the correspondence with the ordinary formalism is achieved by using 

(N\D ai \N) = -±(N\T a <T i \N), (12) 

(Tj being the spin Pauli matrices. 

The scalar form factor in the Skyrme model can be obtained directly from eqs. (§) and ( |l0l) 
and reads 

a Sk {d) = (N\ - C sb {d) \N) = -m 2 Jl [cosF(rf) - 1] . (13) 
On the other hand, the asymptotic form of the chiral angle is determined by C a as pF 
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and hence, for large distances, the leading term in eq. fllBP yields a^(d) = a HB (d). In order 
to test this relationship further, we write 

ml 2 



*L k (d) = mlMN\y/f* -U\N) ^-^tt 
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777 

= -f [(NK\N)(NK\N) + (NK\A)(AK\N)] . (15) 

Using eqs. (|TT1) and (|T^) in the last expression, one concludes that N and A intermediate 
states determine respectively 1/3 and 2/3 of the total value of a^(d). This relative proportion 
is identical to that found recently in the framework of chiral perturbation theory fl25 |, indicating 



that the Skyrme model does provide a rather reliable description of the scalar form factor in 
the heavy baryon limit. 

For systems with B = 2, the standard point of departure for constructing approximate 
solutions is the product ansatz (PA). It uses two undistorted B = 1 hedgehog solutions, whose 
centers are located at two fixed points equidistant from origin along the z axis, so that the 
hedgehog space coordinates are given by y = r + dz/2 and w = r — dz/2. Denoting the 
composite field by U(y,w), one writes 

U(y,w) = U{y)U{w). (16) 

In this configuration, the B = 2 condition is automatically fulfilled, for any distance between 
their centers ||29|| . As the PA keeps the identities of constituent skyrmions, it allows the direct 
incorporation of spin and isospin, through collective rotations of individual hedgehogs. 
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The potential is a function of the distance d and given by 



V(d) 



d 3 r C int (r,dz) 



(17) 



where Ci nt is obtained by using the field U(y,w) in the Skyrme Lagrangian, eqs. and 
subtracting the self energies C[U(y)] and C[U(w)]. This potential works well in the isospin- 
dependent channels, since the OPEP is reproduced for distances larger than 2 fm and it is 
also possible to identify the roles of p and Al mesons |3"0f. On the other hand, problems 



occur in the scalar-isoscalar channel, where the interaction is repulsive at all distances, in sharp 
contradiction with phenomenology. 

Using the definitions F' r = dF(r)/dr, s r = sinF(r) and c r = cosF(r), the central potential 
is given by 
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In order to study its asymptotic structure, we note that the pion fields exist effectively 
only in the neighbourhood of the hedgehog centers. When the distance d is large the skyrmion 
located at (0, 0, d/2) is in the presence of the asymptotic region of U(y), we expand F y , F' y and 
yw around the point w = and write 



\ X X z I X 



(19) 



y ~ 



yw 



—a e 



w. 



X J" 

m^p 2 



x 

2 ^2' 



p 2 + W 2 Z 



W 7 X 



Sm^p 
2x 2 



(20) 
(21) 



where fi, gi are dimensionless polynomials of w z = {z — d/2) and p, which are not displayed 
here. These expressions were tested order by order, by using them in eq. (|T8|) and checking 
that the potential did had the asymptotic structure, as in eq. (^). We found out that it was 
necessary to expand F(y) up to order d~ 2 , in order to have accurate results. 

Replacing eqs. ( |l9l - |2T|) into (18), we obtain an asymptotic contribution of the form 
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for both L a and £4, separately. The values of the parameters K and a, are displayed in table 
0, based on the numerical constants m n = 139 MeV, /„- = 93 MeV and e = 4.0. For the sake of 
comparison, we also present the values of those parameters in the case of the phenomenological 



Argonne potential [31 . 
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U 




Argonne 


K(MeV) 


7.29 


-8.52 


-1.23 


4.80 




0.93 


4.31 


24.4 


1.0 


a 2 


1.78 


12.9 


79.1 


6.0 



Table 2: Coefficients of the multipole expansion of Vc for the product ansatz and the Argonne 
potential, as defined in eq. (|22|) . 

Inspecting this table, one notes that the part of the potential due to L a is attractive, but is 
superseeded by a repulsive contribution coming from the stabilizing term. The net sign of the 
potential is, then, the outcome of a large cancellation. On the other hand, the dependence of 
eq. (|22|) on d is similar to those of both the perturbative chiral calculation, eq. (|6]), and of the 
phenomenological Argonne potential. We stress that this correct geometry is a general feature 
of the model, because it depends only on the form of the B = 1 solution, and not on the specific 
ansatz used to obtain the B = 2 result. In fig. [I] we display the ratios between the full and 
asymptotic PA potentials, as given by eqs. ( |i~8"D and (p2|), with the purpose of illustrating their 
convergence. The interplay between the attractive contribution from L a and the repulsive one 
from £4 can also be seen in fig. 0, where we present the function dVc a (d)/dz, corresponding to 
the integrand in z of expression (|i~8D, for several values of d. One notes that the contributions 
in the neighbourhood of the skyrmion centers are large and positive but, on the other hand, a 
negative region develops as the distance d increases. 

These features of the central potential allow us to identify clearly the stabilizing term as the 
responsible for its repulsive character. Therefore, mechanisms which can reduce the importance 
of £4 may help in producing an attractive interaction. In the next section we discuss a class of 
such mechanisms, associated with deformations of the QCD vacuum. 

3 Constructing B = 2 solutions 

We consider here B = 2 solutions, constructed by using the hedgehog B = 1 skyrmions as 
building blocks, in the framework of the sudden approximation. In general, an ansatz is a 
prescription of the form 

U(y,w) = f[U(y),U(w)], (23) 

where / is a function, chosen according to physical criteria. The construction of such a function 
should follow some guidelines: 

1. the baryon number of the composite configuration must be two for all distances d; 

2. the composite pion field must have the correct quantum numbers, being pseudoscalar, isovec- 
tor and odd under G-parity; 

3. the composite Lagrangian must be chiral symmetric, even under G-parity and invariant 
under the exchange of the two constituent skyrmions. 

The standard constructive approximate solution to the B = 2 system is based on the PA, 
as discussed in the previous section. In this approach, the composite pion field, obtained from 
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0.95 




0.85 



Figure 1: Ratios between the multipole expansion and the exact numerical result for the scalar 
potential (solid) and separate contributions of C a (dashed) and £ 4 (dotted). 



eqs. fllUD and (|16|), is given by 

Ppa = -J- {VyKw + CTw^y ~ ^y x ^w) , (24) 
J TV 

where 7r r is the pion field of the hedgehog with coordinate r and a r = f n cos F(r). The function 

Spa = -J- [VyVw — ^y-^w) ■ (25) 

is the composite analogous of a and satisfies Sp a + P 2 pa = f%. 

The field P pa has a rather serious drawback as a candidate for the pion field, namely 
that it contains an azimuthal term which is both even under G-parity and antisymmetric under 
hedgehog exchange. Hence it does not have good pion quantum numbers, violating requirements 
2 and 3 stated above. 

This motivated us to try to understand whether this problem could be responsible for 
the absence of attraction found in the central potential. We considered several alternative 
possibilities, inspired in the PA. The basic idea is to propose a composite field P, use it to 
define a function S by 

S 2 = fl - P 2 , (26) 

construct the unitary field as 

U=[S + ir-P]/U, (27) 
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Figure 2: Integrand in z of the scalar potential Vc(d), for (a) d — 1 fm, (b) d — 3 fm, (c) 
d = 5 fm, (d) d = 10 fm; vertical axis: arbitrary unity. 



and feed it into the Skyrme Lagrangian. We begin by describing briefly some unsuccessful 
attempts, in order to prevent readers from repeating them. 

The simplest exchange-symmetric ansatz would be the average P = (7v y + ir w ) /2. However, 
when d = 0, one has F y = F w = F r and hence P = f n sin F r r corresponds to a B = 1 field, 
which must be disregarded. 

This suggests that, in order to obtain B — 2, it is mandatory to mix 7r and a. In the case 
of the PA, which yields B = 2 at all distances, we note that the chiral constraint between P pa 
and S pa allows one to write 

U pa = [S pa + ir-P pa ] IU = e- lT - UF ^ , (28) 

where u is a unit vector, taken as pointing always away from the origin of the coordinate 
system, and F pa is a profile function. In fig. |3] we display the behaviour of this angle along 
the axes z and p, for various values of the internucleon distance d. The solid lines correspond 
to the case d = 0, which is spherically symmetric and it is possible to see that, along both 
directions, the chiral angle varies smoothly from 2n at the origin to at infinity. In the case 
d = 0.5 fm, shown in dotted lines, one notes that a discontinuity has appeared along the p axis. 
This discontinuity increases with distance and, at d cr u = 0.86 fm, the chiral angle is such that 
F pa {.P = 0, z — > 0) = 2n and F pa (p — > 0, z = 0) = 0. Therefore, at this critical point, it is more 
natural to set F pa (0, 0) = and to work with two separate solutions, such as illustrated by the 
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Figure 3: Profile function F for the product ansatz, in unities of tt, along the z (left) and p 
(right) axes, for various values of d. 

dashed and dot-dashed lines, corresponding to 0.9 and 2.0 fm, respectively. This suggests that, 
from d cr it onwards, each of the interacting skyrmions acquires a considerable individuality. 
The combination 

P = — {vyitw + o w -K y ) (29) 

is interesting, for it has an explicit physical meaning. As the function cr(r) is associated with 
the quark condensate that surrounds the baryon labeled by r, this field P represents each 
skyrmion immersed in the distorted vacuum of the other one. The condition allows one 
to determine S up to a sign. In the case B = 1 the field a changes sign when one goes from 
infinity to the origin and the same happens when B = 2. The sign of S is also important and, 
in order to fix it, we note that the behaviours of eqs. (p4|) and ( f29|) along the z axis are identical, 
since the azimuthal component vanishes. We then forced the condition S = S pa along this axis. 
However, this ansatz, based on eq. ( p9|) , gives rise to a baryon number which varies with d, as 
shown in fig. f|, and had to be abandoned. 

This discussion illustrates the fact that it is not trivial to build an ansatz with a good 
topology. We thus decided to adopt simultaneously the pion field as given by eq. fl29"D and the 
function S ap of the PA, eq. Q2"5|), for its topology is automatically correct. With this option, 
the unitarity constraint reads 

Sl P + P 2 = flv 2 , (30) 

where 
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Figure 4: Baryon number as a function of distance d (fm), for the ansatz given by eqs. 
and (M). 



and the pion field becomes in fact P/r]. This form for the dynamical variable is the same as 
that proposed by Nyman and Riska, in their symmetrized product ansatz [|17]]. This ansatz has 
a topology similar to the PA, as illustrated in fig. [j| The corresponding baryon number density 
is given in appendix A and, in the classical case, yields B = 2 for all distances when integrated 
over space, as shown in fig. ||. 

4 Results and conclusions 



In order to derive the potential, we use the quantized fields of eq. obtained by rotat- 

ing the constituent skyrmions. This idea of rotating individual hedgehogs corresponds to an 
approximation and deserves some attention. The quantization of a hedgehog, as discussed in 
sect. ^ amounts to multiplying the classical field by the matrix D, which depends on three 
free parameters. In the case B = 1, this procedure does not change the baryon current. This 
can be seen by writing the baryon density for quantized fields as 



and using the result 



in order to obtain 



B° — — €abc €a f 3 "t daDaiTTi O^DpjTXj d c D^ k TT k 
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(33) 
(34) 
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Figure 5: Baryon number as a function of the separation distance (fm) for the symmetrized 
ansatz (solid); separate contributions B 1 (dashed) and B 2 (dotted), as defined in eq. |39|. 

This shows that the baryon density is the same for both quantized and classical fields. 

Analogously, in the case B = 2, quantization would require a matrix D, depending on 
six collective coordinates. However, the determination of this general matrix may prove to 
be very difficult and, in the spirit of the the sudden approximation, one normally uses D pa 
E>(y)lM + i(y)£)M : wne re I is an identity matrix and £>M are operators over the skyrmions 
labeled by y and w respectively. The price one pays for this approximation is that it leads to a 
quantized baryon number which depends on and D^ w \ This happens because the relation 
equivalent to eq. fl33|) does not hold for the approximate matrix D and hence does not represent 
a major shortcoming for the symmetrized ansatz (SA). Indeed, as pointed out by Sternheim 
and Kalbermann jTjl, this poses problems for short distances only. 

A collective rotation of the pion field P, as in the B = 1 case, would leave S/rj unmodified, 
as a classical function. However, this would also mean to treat the scalar product iZy ■ tc w as a 
classical quantity and would lead to serious contradictions, for the OPEP content of the isospin 
dependent channels relies on the quantum character of such a scalar product in eq. (|25|) . 

Therefore, the individual rotation of each pion field is more consistent with a constructive 
approach, although not free of problems. In principle, every pion field 7r in the composite 
Skyrme Lagrangian should be quantized. When applying this prescription to the dynamical 
variable of the SA, one has to deal with the functions t]~ 2 and r? -4 , which depend on pion fields, 
coupled to operators D. The meaning of the quantized r?™ 2 is that of a power series in D, which 
involves products of arbitrary numbers of these matrices and hence can only be handled by 
resorting to truncation. With this limitation in mind, we treat rj~ 2 as a polynomial in 7r and 
thus its expectation value between two-nucleon states can be evaluated without ambiguities. 
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In order to test the implications of this assumption, in the sequence we present results with 
two versions of r/, namely, a classical one, 

v7 2 = {i - [i - [y-™?] 4*1 }~\ ( 35 ) 

and a quantized one, truncated at the first order in the D expansion, given by 
(NN\rjf\NN) = (NN\{l- S y w + D^D^D$D$ M j (* v ) k (tt w ), /£ \NN) 
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Replacing the pion field of the symmetrized ansatz into the interaction Lagrangian used to 
calculate the potential, one has 

C , — _i_ n(y) nO) /o 7 \ 

where the labels S and V stand respectively for isoscalar and isospin dependent parts of Ci nt . 
Using this result in eq. fll7|), one gets 



V{d) = V C + t^-t (w) [trM-trWVss + S 12 V T } , (38) 

where Vc, Vss and Vp are the usual central, spin-spin and tensor components. All terms receive 
both G-parity odd and even contributions. 

The G-parity odd components of spin-spin and tensor terms of the potential are shown in 
fig. H, for the two possible choices of r], compared to the PA and pure OPEP results. One sees 
that all curves coincide for distances larger than 2 fm, indicating that all ansatze reproduce 
asymptotically the OPEP. The results for the G-parity even terms are of minor importance 
here, as we are interested in the long range behaviour of the potential, but they are included 
for completeness in fig. 0. One should note that in the case of rj c , the potentials present a 
singularity at d ~ 0.6 fm, due to a root of r) c . 

Results for the scalar component Vc are presented in figs. |8| and ^. In the former we display 
the behaviour of the PA (left) and the predictions from the SA with r\ — 1 (right), which is 
non-unitary and considered just for pedagogical purposes. Inspecting it one learns that the SA 
includes a contribution from £ 2 , that was not present in the PA. Moreover the contribution 
from £ 4 is negative, and so is the net result for Vc- 

The two valid options for the SA considered here, based on t] c and T] q , are given in fig. |9|. 
In both cases we observe that the unitarity constraint restores the repulsion due to £4, but in 
such a way that the net result is asymptotically attractive. On the other hand, the amount of 
overall attraction found in the SA depends on the specific quantization prescription adopted. 
At very large distances, the curves corresponding to rj c and r\ q have the same geometry and 
yield respectively the following approximate values for the intensity of the potential: K c = 14 
MeV and K q = 57 MeV. Comparing them with the empirical values in table [1], it is possible to 
see that predictions from the SA are qualitatively reasonable. 

In summary, we have shown that the SA provides the correct baryon number for the two- 
nucleon system in the Skyrme model, as well as attractive central iViV potentials. The correct 
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Figure 6: Comparison among Skyrme model G-parity odd spin-spin (left) and tensor (right) 
potentials from symmetrized ansatz, with classical (dotted) and quantized (dashed) versions of 
T), from product ansatz (long dashed) and OPEP (solid). 
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Figure 7: G-parity even spin-spin (left) and tensor (right) potentials from symmetrized ansatz, 
with classical (dotted) and quantized (dashed) versions of 77, and from product ansatz (long 
dashed). 

quantitative feature is somewhere between the values obtained for the two versions of the 
normalization function rj. 

We conclude that it is indeed relevant to the central potential to eliminate the term with a 
wrong G-parity from P pa . We expect that a deeper and more careful study of the quantization 
procedure will lead to a more accurate evaluation of the amount of attraction coming from the 
SA. 
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A Appendix: The baryon number 

The explicit calculation of the zero-component of the baryon current yields 



B Q = B° X +B Q 21 



where 



B° 2 (r;d)- 



2tt 2 (c y c w - (y-w) s y s w ) r] 



2 °» 



J w . a y 



W 



I - (y-w) 



^y^w I ^y-L y 



c„c m [ c„H - °JL \{ CwF ' w - ^ ) - slslFl.F. 



w 



.2 n 2 tt' I?' 



y w y w 



Sy$w 



J y ( [ c v~^ + Cw ~^) (° wSy + (y'™) c v Sw ) + ( l ~ (y ' ^) 2 ) s y°y ( F w - c 



w 



(39) 



[y <-> w) 



with J y = ((£r • wf - l) c y F y s w + (y-w) (s y ^- -{y-w) s w s f). 

The numerical integration of B® is tricky due to the presence of the function S in the 
denominator. Results are shown in fig. |5|, as functions of separation distance d. It shows that 
the SA presents the correct topology for the iViV system. 



References 

[1] W. N. Cottingham and R. Vinh Mau, Phys. Rev.lZQ, 735 (1963); 

W. N. Cottingham, M. Lacombe, B. Loiseau, J. M. Richard and R. Vinh Mau, Phys. Rev. 
D 8, 800 (1973); 

M. Lacombe, B. Loiseau, J. M. Richard, R. Vinh Mau, J. Cote, P. Pires and R. de Tourreil, 
Phys. Rev. C 21, 861 (1980). 

[2] R. Machleidt, K. Holinde and C. Elster, Phys. Rep. 149, 1 (1987). 

[3] G. E. Brown and J. W. Durso, Phys. Lett. B 35, 120 (1971). 

[4] M. Chemtob, J. W. Durso and D. O. Riska, Nucl. Phys. B38, 141 (1972). 

[5] M. Gell-Mann and M. Levy, N. Cim. 16, 705 (1960). 

[6] T. H. R. Skyrme, Nucl. Phys. 31, 556 (1962). 

[7] E. Witten, Nucl. Phys. B160, 57 (1979). 

[8] G. S. Adkins, C. R. Nappi and E. Witten, Nucl. Phys. B228, 552 (1983); 
G. S. Adkins and C. R. Nappi, Nucl Phys. B233, 109 (1984). 



16 



[9] S. Weinberg, Phys. Lett. B 251, 288 (1990); Nucl. Phys. B363, 3 (1991). 

[10] C. Ordonez and U. van Kolck, Phys. Lett. B 291, 459 (1992); 

L. S. Celenza, A. Pantziris and C. M. Shakin, Phys. Rev. C 46, 2213 (1992); 

J. L. Friar and S. A. Coon, Phys. Rev. C 49, 1272 (1994); 

C. A. da Rocha and M. R. Robilotta, Phys. Rev. C 49, 1818 (1994); 

M. C. Birse, Phys. Rev. C49, 2212 (1994). 

[11] G. Holder, in Numerical Data and Functional Relationships in Science and Technology, 
edited by H. Schopper, Landolt-Bernstein, New Series, Group I, Vol. 9b, Part 2 (Springer- 
Verlag, Berlin, 1983). 

[12] C. Ordonez, L. Ray and U. van Kolck, Phys. Rev. Lett. 72, 1982 (1994); Phys. Rev. C 53, 
2086 (1996). 

[13] M. R. Robilotta and C. A. da Rocha, Nucl. Phys. A615, 391 (1997); 

J.-L. Ballot, C. A. da Rocha and M. R. Robilotta, Phys. Rev. C 57, 1574 (1998). 

[14] N. Kaiser, R. Brockman and W. Weise, Nucl. Phys. A625, 758 (1997); 
N. Kaiser, S. Gerstendorfer and W. Weise, Nucl. Phys. A637, 395 (1998). 

[15] A. Jackson, A. D. Jackson and V. Pasquier, Nucl. Phys. A432, 567 (1985). 

[16] R. Vinh Mau, M. Lacombe, B. Loiseau, W. N. Cottingham and P. Lisboa, Phys. Lett. B 
150 259. (1985). 

[17] M. Nyman and D. O. Riska, Phys. Lett. B 203, 13 (1988). 

[18] M. Sternheim and G. Kalbermann, Phys. Lett. B 225, 27 (1989). 

[19] T. S. Walhout and J. Wambach, Int. J. Mod. Phys. E 1, 665 (1992). 

[20] J. J. M. Verbaarschot, T. S. Walhout, J. Wambach and H. W. Wyld, Nucl. Phys. A461, 
603 (1986); 

V. B. Kopeliovich and B. E. Shtern, JETP Lett. 45, 203 (1987). 
[21] T. S. Walhout and J. Wambach, Phys. Rev. Lett. 67, 314 (1991). 
[22] J. M. Eisenberg and G. Kalbermann, Int. J. Mod. Phys. E 5, 423 (1996). 
[23] T. Gisiger and M. B. Paranjape, Phys. Rep. 306, 109 (1998). 

[24] C. M. Maekawa, J. C. Pupin and M. R. Robilotta, Phys. Rev. C 61, 064002 (2000). 
[25] M. R. Robilotta, submitted for publication. 

[26] W. B. Kaufmann and G. E. Hite, Phys. Rev. C 60, 055204 (1999). 



17 



[27] M. M. Pavan, R. A. Arndt, I. I. Strakowsky and R. L. Workman, nN Newsletter, 15, 118 
(1999). 

[28] G. S. Adkins, in Chiral Solitons, ed. F. Liu, World Scientific, p. 99 (1987). 

[29] I. Zahed and G. E. Brown, Phys. Rep. 142, 1 (1986). 

[30] H. Yabu and K. Andou, Progr. Theor. Phys. 74, 750 (1985). 

[31] R. B. Wiringa, R. A. Smith and T. L. Ainsworth, Phys. Rev. C 29, 1207 (1984). 



18 



